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We study a system of D-branes localized near an intersection of Neveu-Schwarz fivebranes, 
that is known to exhibit a landscape of supersymmetric and (metastable) supersymmetry 
breaking vacua. We show that early universe cosmology drives it to a particular long-lived 
supersymmetry breaking ground state. 



1. Introduction 

String theory is believed to give rise to a rich landscape of supersymmetric and non- 
supersymmetric vacua. A natural question is whether early universe cosmology can provide 
a dynamical selection mechanism among these vacua. The answer is unknown in general, 
both because the full structure of the landscape is at present not well understood, and 
since real time dynamics in string theory is often hard to study. 

To gain insight into this issue, it might be useful to consider simplified versions of the 
landscape, which are expected to play a role in the full problem, but are more amenable 
to a detailed analysis. The purpose of this paper is to discuss an example of a "mini 
landscape" of this sort, focusing on the question of dynamical vacuum selection. 

We will consider a system of Neveu-Schwarz fivebranes connected by Dirichlet four- 
branes, which is known to be useful for realizing supersymmetric gauge theories in string 
theory and has been used in the last few years to study spontaneous supersymmetry 
breaking [0||]. In a particular region of its parameter space, the low energy dynamics 
of this brane system reduces to a Wess-Zumino (WZ) model which was analyzed in 0. 
In another, one has to take into account classical string effects, which give rise to non- 
renormalizable terms in the low energy effective Lagrangian, suppressed by a relatively 
low UV scale [0,0. We will refer to the two regions as the field theory and string theory 
regimes, respectively. 

Interestingly, the vacuum structure found in || in the field theory regime is very similar 
to that obtained in the string theory one in 0,0]. In both cases one finds a metastable 
supersymmetry breaking ground state, located at the origin of an approximate moduli 
space, and separated by a wide potential barrier from a supersymmetric ground state. 
The effects that lift the moduli space are different in the two regions. In the field theory 
regime they are due to a one loop contribution to the effective potential of the low energy 
WZ model, while in the string theory regime their origin is the classical gravitational 
interaction between D-branes and iVS'S-branes. 

It is natural to ask whether early universe dynamics drives the system to the true 
(supersymmetric) vacuum, or to the metastable (supersymmetry breaking) one. In the field 
theory regime this question was addressed in . These authors analyzed the structure 
of the finite temperature effective potential and argued that, under mild assumptions, the 
system ends up in the metastable vacuum. 

In this paper we will study this question in the string theory regime. We will see 
that early universe dynamics drives the system to the metastable state, like in the field 



1 



theory regime, but the details are different. Furthermore, in string theory one encounters 
an important issue that does not arise in the field theory discussion. In studying the 
early universe dynamics in field theory, the couplings in the Lagrangian are taken to be 
fixed (i.e. time independent ).0 In string theory, these couplings are parameters of the 
brane configuration, and may depend on time. We find that the time evolution of these 
parameters leads in general to a singular limit of the ISS brane configuration. 

Our main purpose in this paper is to study the generalized ISS model which, as we 
will see, does not suffer from the above problem. The vacuum structure of this model 



was analyzed in |T2| in the string theory regime and in |T^,|n| in the field theory one. It 
explicitly breaks R-symmetry, and thus is promising for phenomenological applications. 
The field theory analysis leads in this case to a number of supersymmetric and metastable 
non-supersymmetric vacua, which are generalizations of those found in ||. In the string 
regime one finds many additional metastable vacua, which are qualitatively different from 
those of H and potentially interesting for phenomenology. We will see that early universe 
dynamics provides a selection mechanism among the stable and metastable vacua of the 
generalized ISS brane system. 

The plan of this paper is the following. We start in section 2 by reviewing some field 
theoretic results. We describe the ISS model, its R-symmetry breaking generalization, and 



the results of [p|-p!T| on the early universe evolution of this model. In section 3 we review the 
brane construction of the (generalized) ISS model and the role of the various parameters 
that appear in it. In section 4 we study the dynamics of the brane systems of section 3 in 
typical excited states that are expected to play a role in the early universe. We argue that 
the main effect of the excitation is to turn the Neveu-Schwarz fivebranes into non-extremal 
black branes with non-zero horizon radii. The light translational modes on the D-branes 
are excited as well. Thus, the study of early universe dynamics of the brane systems in 
the classical string theory regime amounts to analyzing the dynamics of the .D-branes in 
the background of non-extremal fivebranes. 

Applying these considerations to the ISS brane system, we find that at early times 
the system is in a phase with unbroken gauge symmetry. As the fivebranes relax to their 
extremal limit, the system undergoes a phase transition to the metastable state. This 



conclusion is similar to that found in the field theory regime in [p HTlJI , but the details are 
different. We also point out that the attraction of the non-extremal AS'S-branes to each 



1 Except for the usual RG evolution with the scale. 
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other naturally leads the ISS system at late times to a state in which some of the branes 
coalesce. Hence the system of is an unnatural endpoint of this dynamics. We then 
turn to the generalized ISS brane configuration which, as mentioned above, does not suffer 
from this problem. The rich vacuum structure of this system makes the problem of vacuum 
selection particularly interesting. We show that the dynamics of the D-branes in the 
background of non-extremal fivebranes drives the system to a specific non-supersymmetric 
ground state. 

In section 5 we comment on our results and possible extensions. Two appendices 
contain the details of some of the calculations. 



2. Review of field theory results 

In this section, we briefly review the ISS model its finite temperature extension 
P||TT|], and the generalization discussed in flT3|jnj] . The main purpose of this discussion is 



to set the stage for the subsequent analysis. 



2.1. The ISS model 



Consider Af = 1 supersymmetric QCD (SQCD) with gauge group SU(N C ) and Nf 
flavors of (anti) fundamentals Qi, Q l , % = 1, 2, • • • , Nf. This ("electric") theory is believed 
to be Seiberg dual [ffBfl to a "magnetic" gauge theory, with gauge group SU (Nf — N c ), Nf 



fundamentals q l , qi, and a gauge singlet meson field which is coupled to the magnetic 
quarks q, q via the superpotential W m = hq&q. 

The two theories are expected to be equivalent in the infrared, i.e. well below their 
dynamically generated scales. When one of them is weakly coupled at long distances, 
the other is strongly coupled. In the range N c < Nf < 3N c /2, where the analysis of 
takes place, the electric theory is strongly coupled in the infrared, while the magnetic 
one is infrared free. In this regime, the electric theory can be viewed as providing a UV 
completion of the non-asymptotically free magnetic theory. 

To study spontaneous supersymmetry breaking, one needs to turn on a small mass 
for the electric quarks Q, Q. In the magnetic theory, this corresponds to deforming the 
superpotential as follows: 

W m = hq i &q i -hn 2 & i . (2.1) 



The superfields $, q, q are normalized such that their Kahler potential is canonical, up to 
corrections suppressed by powers of the dynamically generated scales. Thus, the classical 
bosonic potential takes the form 

V = \h\ 2 (\qq-(j?I Nf \ 2 + \q$\ 2 + \$q\ 2 ) , (2.2) 

where 1^. is the Nf x Nf identity matrix. As in the O'Raifeartaigh model, the potential 
( |2.2| ) cannot be set to zero due to the rank of the matrix qq; hence supersymmetry is 
broken. The minimum of the potential corresponds (up to global symmetries) to 
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The N c x N c matrix valued field X has a flat potential, and parameterizes a pseudo-moduli 
space of classical vacua. Since supersymmetry is broken, the fate of this moduli space is 
sensitive to small corrections to the dynamics. In the regime of interest for the analysis 
of 0, the leading correction is the one loop Coleman- Weinberg potential, which behaves 
near the origin like 

Vi = l ^r^{N f -N c )\h 2 ^\ 2 TT^X + -- - . (2.4) 

It gives a mass of order \h 2 n\ to the pseudo-moduli, and stabilizes them at the origin. Non- 
perturbative gauge dynamics gives rise to additional, supersymmetric, vacua at q = q = 
and large $, 



<*>=*°=fUJ ,N " (2 - 5) 

where A is the dynamically generated scale of the magnetic gauge theory. The super- 
symmetric vacua ( [2.5| ) are separated from the non-supersymmetric one at the origin by a 
potential barrier. The vacuum at the origin is long-lived for \i -C A. 

2.2. Finite temperature 



The authors of [§|-|TT1 asked whether early universe dynamics is likely to drive the ISS 
theory to the metastable state at $ = 0, or to the lower energy supersymmetric state 
(|2.5|) . They assumed that after inflation and reheating the universe is in a thermal state 
with an adiabatically decreasing temperature.ll In that case, one needs to evaluate the 



2 For this to be true, the thermal equilibration time must be much shorter than other relevant 
time scales, such as the rate of change of the temperature, and the Hubble time. 
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finite temperature corrections to the effective potential Vq + V±, (|2.2| ), (2~1). The thermal 
effective potential Vr, which depends non-trivially on both the pseudo-moduli X, (|2.3|), 
and magnetic quarks, g, g, takes the form 

V T = V + V 1 + VP, (2.6) 

where [1TB 



VP (q, q,X) = ^J2 ±n * d PP 2 ln i 1 ^ exp(- ^p 2 + m 2 (g, g, A)/T^ . (2.7) 

Here rij and m; label the degeneracies and masses of different fields, respectively. The sum 
is over all the fields in the theory, with the upper sign for bosons and the lower one for 
fermions. 

For temperatures slightly above 

/( (2.8) 



' V N c + 2N f ' 

the effective potential ( [2 .6] ) has a global minimum at the origin of field space, q = q = 
$ = 0, and a local one at the location of the supersymmetric vacuum, ( |2.5| ). The global 
and local minima are separated by a potential barrier, just like in the zero temperature 
case 0. At these temperatures, thermal fluctuations are not sufficient to push the system 
over the barrier. Let us assume for now that at that point the system is in the basin of 
attraction of the global minimum at the origin (we will revisit this assumption below), and 
consider its evolution. 

When the temperature decreases below T c , (|2.8|) , the origin becomes unstable to con- 
densation of the magnetic quarks g, g, and the system undergoes a second order phase 
transition to a phase with broken gauge symmetry, which for low temperatures approaches 
the metastable vacuum ( |2.3| ). Initially, the supersymmetric vacuum still has a higher free 
energy than the metastable one, but when the temperature drops below a value of or- 
der \/hT c (we are assuming that the Yukawa coupling h 1), the free energy of the 
supersymmetric vacuum becomes lower. 

At that point, the metastable vacuum can decay to the supersymmetric one by tun- 
neling. However, the time scale for this first order phase transition is governed by the 
parameter yu/A, which is taken to be small to ensure the metastability of the vacuum ( |2.3| ) 
at zero temperature. This ensures that the metastable vacuum does not decay in the early 
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universe as well, at least when the temperature is of order T c ( |2.8|) or lower. Thus, at late 
times, when the temperature T —>■ 0, the system remains in the metastable vacuum. 

In the above discussion we assumed that at T ~ T c the system is near the metastable 
vacuum at the origin. If, instead, it is near the super symmetric vacuum, it would remain 
there at late times. Thus, an important question is whether dynamics in the regime T > T c 
leads predominantly to one or the other outcome. 

The authors of |J argued that high temperature dynamics drives the system to the 
metastable vacuum. The basic point is that as the temperature increases, the supersym- 
metric vacuum becomes less stable, until it disappears at a temperature of orderi h&o, 
( [2 . 5| ) . If the reheating temperature after inflation is higher than this value, $ rolls towards 
the origin and is likely to get trapped there as the temperature decreases (as described 
above) . 

We see that the ISS system provides an example of dynamical vacuum selection. At 
zero temperature it has two (types of) ground states, one of which is stable and super- 
symmetric and the other metastable and non-supersymmetric. One might expect that the 
system will end up at late times in the supersymmetric ground state, or assign some other 
probabilities to the two vacua. Instead, one finds that dynamics favors the metastable vac- 
uum over the supersymmetric one. Below, we will generalize this picture to brane systems 
in string theory, which, as mentioned above, have in some cases a more intricate vacuum 
structure. 



2. 3. Generalized ISS model 

The unbroken R-symmetry of the ISS model leads to anomalously small gaugino 



masses in the context of gauge mediation (see e.g. ]l7j for a recent discussion). Thus, 
it is interesting to consider R-symmetry breaking deformations of this model. An example 
of such a deformation, studied in []T3|JT^1 (see also [p~8| , pT9f| ) , is obtained by replacing the 
magnetic superpotential ( |2.1|) by 

1 



W m = hq$q- hfi 2 Tr$ + -h 2 ^Tr$ 2 



which gives rise to the bosonic potential 



V = N 2 {\qq~ (J?I Nf + h^\ 2 + \q0>\ 2 + \$q\ 2 ) 



(2.9) 



(2.10) 



3 In Jllj j it was suggested that the temperature at which the SUSY vacuum disappears is 
somewhat lower than h&o, but still much higher than T c , fl2.8|) . 
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We will restrict below to the case <C fj,, for which one can consider the last term in 
( |2.9|) as a small perturbation of $Z7\\). 

To study the vacuum structure of the theory with deformed superpotential ( |2.9| ), it is 
convenient to parameterize the light fields as follows: 

/«* \ (n 2 h o \ 

h$ = \ h$ n I o ^ . (2.11) 

\0 £l Nf _ k _ n J \ Nf - k - n J 

Here 99 and are n x (iVj — N c — k) dimensional matrices, while $ n and are n x n 
matrices. This parameterization is motivated by the brane construction, which will be 
reviewed in the next section. 

Unlike the case = 0, the superpotential ( |2.9| ) admits classical supersymmetric 
vacua,@ which correspond in (|2.11| ) to n = 0, k = 0, 1, • • • , Nf — N c . In these vacua, the 
SU(Nf — N c ) magnetic gauge symmetry is broken to SU(Nf — N c — k), by the expectation 
value of the magnetic quarks. Some components of $ get a large (3> //) expectation value. 
The fields $ n , 99, <p are absent since n = 0. 

There are also metastable vacua, in sectors with n^O. In studying them, it is useful 
to discuss separately two cases: 

(1) For k = Nf — N c the magnetic gauge group is completely broken. The fields 99, <p 
are absent in this case, so the n x n matrix tptp in (|2.11|) should be set to zero. The 
potential for $ n has two contributions. One is the classical potential ( ^.10 ), which 
pushes it towards the supersymmetric vacuum /i$ n = (// 2 / ' \i^)l n . The other comes 



from the one loop potential (|2.4|) , and provides an attractive force towards the origin. 



Balancing the two against each other, one finds |13(] a local minimum of the potential 
near the origin, at 



In order for this vacuum to be in the regime of validity of the approximation ( |2.4j ), 
must be very small, hfi(Nf — N c ). 

(2) For k < Nf — N Cl part of the SU(Nf — N c ) gauge group is unbroken by ( |2.11| ), 
and the quark fields 99, <p need to be included in the discussion. Looking back at 
( |2.10|) one finds that they are tachyonic, unless the minimum of the potential for the 



For some values of Nf, N c , there are additional supersymmetric vacua that are quantum in 
nature [13]. 
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pseudo-moduli Q2.12 ) is at $ n > \i. In this regime, the approximation to the one loop 
potential ( |2.4j) is not valid. The full one loop potential has no metastable vacua in 
these sectors |14| (see also [^8j). 
As shown in [[12] , the brane system that realizes the generalized ISS field theory in a certain 
regime in its parameter space has a similar phase structure, but it has additional metastable 
states in sectors with k < Nf — N c . In the next section we review this construction and 
its vacuum structure. 



The discussion of [p|-|TT|| of finite temperature dynamics of the ISS model is easy to 
generalize to the case of non-zero \x$. The vicinity of the metastable vacuum ( |2.12| ) can be 
studied as in the previous subsection, since the small \i$ does not modify the potential by 



much (see |20[ for a related discussion). The fate of the supersymmetric vacua discussed 
above, (|2 .11 ) with n = 0, can also be studied. The situation here is better than in ||, since 
these vacua are visible classically in the WZ model, and do not require non-perturbative 
effects in the non-asymptotically free magnetic theory. By evaluating the effective potential 
fl2.7p , one finds that the SUSY vacuum disappears when the temperature exceeds a value 
of order /j, 2 / [i^ (with a small coefficient, of order 0.01 — 0.1 for the values of the parameters 
that we considered). 



3. Review of the brane construction 

In this section we review the brane construction of the ISS model P-EL and its general- 



ization presented in [12]. To describe it, it is convenient to decompose the 9 + 1 dimensional 



spacetime of type IIA string theory as follows: 

1R 9 ' 1 = IR 3 ' 1 xC„ x <C W x M y x JR X 7 , (3.1) 

with 

v = x A + ix 5 , w = x 8 + ix 9 , y = x 6 . (3.2) 

The brane configuration we consider is depicted in figure 1. All branes are extended in 
the IR 3 ' 1 labeled by ( ). The iVS'S-branes denoted by NS and NS' are further 

extended in v and w, respectively. The L>4-branes are stretched between other branes as 
indicated in the figure. 
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Fig. 1: The ISS brane configuration. 

The fivebranes alone preserve N = 2 supersymmetry in four dimensions. They can be 
thought of as providing a dual description of a non-compact Calabi-Yau (CY) geometry. 
For yi = yisrs = v % = this geometry is singular, and can be described as in |2TJ. Varying 
these parameters corresponds to smoothing out the singularity. The D4-branes stretched 
between iV,S5-branes in figure 1 correspond in the dual picture to L>5-branes wrapping 
small two-cycles in the non-compact CY. The CY description is useful when one of the 
directions transverse to the fivebranes is compactified on a small circle. We will consider 
the opposite limit, in which the fivebrane picture is the appropriate one. 

At low energies, the brane configuration of figure 1 reduces to an N = 1 supersym- 
metric gauge theory with gauge group SU(Nf — N c ) x SU(Nf) x U(l), Nf flavors of 
bifundamentals q, q with opposite U(l) charges, and an adjoint of SU(Nf), $. In the 
brane picture, the eigenvalues of $ correspond to the locations in the u>-plane of the D4- 
branes stretched between iV£"-branes. The precise relation between $ and w is described 
e.g. in 0. 

The chiral superfields q, q and $ have canonical KaMer potential, and superpotential 



This theory is essentially the same as that described in section 2.1, except for the gauging 
of the flavor symmetry group. The theory in which the flavor group is not gauged can also 
be realized as a low energy theory on branes, by replacing the iVS^-brane in figure 1 by 



mag 



hq&q - /i/i 2 Tr$ . 



(3.3) 
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Nf D6-branes [22]. For our purposes below it will be convenient to consider the system 
with gauged flavor group. 

The parameters in the low energy theory on the branes can be expressed in terms of 
string theory quantities as follows. The gauge couplings g\ M for the different factors of 
the gauge group are proportional to g 8 divided by the extent in y of the relevant fourbranes 
0. The parameters in the superpotential ( [OP are given by 



h 2 



S7r 2 g s l s 
Ay 



V2 



(3.4) 



where Ay = y% — yi, and we made a choice of phase for the chiral superfields such that 
the complex couplings h and fx are real and positive. 

The geometric description is reliable when the distances between the various branes 
are large (relative to Z s ), and the string coupling g s is small. In this regime, the gauge 
couplings and Yukawa coupling h are small. The mass [i, which is determined by the 
geometric parameter i>2, is typically above the string scale. If \x is sufficiently small, the 
low energy dynamics of the branes is well described by magnetic SQCD. In general, the 
low energy effective Lagrangian receives contributions from other sources [^,0. 
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Fig. 2: The marginally stable brane configuration. 



The brane configuration of figure 1 is unstable to reconnection of the Nf — N c "color" 
D4-branes with some of the Nf "flavor" D4-branes, leading to that of figure 2. The 
resulting configuration is marginally stable. It contains an iV c x N c matrix X of massless 
fields with a flat potential. The (pseudo-) moduli space labeled by X is a direct analog of 
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that found in the low energy field theory analysis in section (2.1). As discussed after eq. 
( p.3|) , the exact flatness of the potential for X may be spoiled by small corrections to the 
dynamics. 

In the field theory regime {i.e. when the geometric parameter v% in figure 2 is small, 
so that fi <C m s ), the leading correction is the CW potential ( |2.4j) , which lifts the moduli 
space by giving a mass to X. In string theory language, the CW potential (|2.4j ) is a 
string one loop effect, and is thus suppressed in the g s — > limit. In the regime where all 
distances in figure 2 are large in string units, the dominant contribution to the potential for 
X comes from the classical gravitational attraction of the iV c L>4-branes to the A^-brane 
||. From the point of view of the low energy effective field theory for X, this attraction 
gives rise to non-renormalizable D-terms that, together with the linear superpotential 
( p.l|) , provide a potential which stabilizes X at the origin. 

Thus, the fate of the pseudo-moduli space labeled by X is the same in the field theory 
and string theory regimes - it is replaced by an isolated supersymmetry breaking vacuum 
at the origin - but the mechanisms are different. The supersymmetric vacuum ( [2.5| ) is not 
visible in the classical brane description. In string theory it owes its existence to quantum 
(g s ) effects, and is thus located at a very large value of w. 

So far in this section we reviewed the realization of the ISS model as a low energy 
theory on branes. Since our main interest is in the generalized ISS model described in 
section 2.3, we next turn to the corresponding brane system. The brane configuration, 
which is depicted in figure 3, is obtained from that of figure 2 by a rotation of the AS^- 
brane by an angle 9 in (v, w). The parameters h, \x are again given by (|3.4| ), while 



tan# 

The orientation of the AS^-brane leads to an important difference in the dynamics of the 
configuration of figure 3 compared to that of figure 2. While the latter has (to leading 
order) a pseudo-moduli space of supersymmetry breaking vacua labeled by the positions 
of the N c D4-branes in w, in the former this moduli space is lifted and replaced by a 
supersymmetric vacuum, in which the L>4-branes are located at the point (v, w) = (0, wo), 
with 

wq = i>2 cot 9 . (3-6) 

At that point, the projections of the A^S'-branes on the (v, w) plane intersect, and the 
D4-branes stretched between them preserve N = 1 supersymmetry. 
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Fig. 3: The deformed ISS brane configuration. 
In fact, the brane system of figure 3 has many supersymmetric vacua, which were 



described and compared to field theory in ||12|| . From the perspective of the WZ model 
with superpotential (|2.9|) , these vacua are obtained by analyzing the zeroes of the bosonic 



potential (|2.10[) . The brane system also has many non-supersymmetric vacua [[12 



The classical supersymmetric and non-supersymmetric vacua are exhibited in figure 
4. They are labeled by two integers k, n, which are the same as those defined in the 
field theory in eq. (|2.11| ). The Nf — N c — k D4-branes stretched between the NS and 
iV5^-branes give rise to an unbroken SU(Nf — N c — k) subgroup of the magnetic gauge 
group. A SU(k) x SU(n) x SU(Nf — k — n) subgroup of the gauge group is unbroken as 
well. 

In the low energy field theory, these vacua can be understood in terms of the decom- 



position (|2.11| ). As mentioned there, the supersymmetric vacua correspond to n = 0. They 
can be obtained by moving the n L>4-branes in figure 4 to the point (v, w) = (0, wq), where 
they join the Nf — k — n D4-branes already there. 

The vacua with n > seem at first sight unstable, since the n D4-branes can decrease 
their energy by moving towards the aforementioned supersymmetric vacua. However, the 
attraction of the D4-branes towards the iVS'-brane provides a counterbalancing force. The 
combination of the two effects leads to the existence of a local minimum of the energy 



at a value of w which was computed in ||12|| , and interpreted from the low energy field 
theory point of view in @ . These non-supersymmetric vacua can be made parametrically 
long-lived by tuning g s and other parameters of the brane configuration. 
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Fig. 4: The vacuum structure of the deformed brane configuration. 

As is clear from figure 4, there are two qualitatively different types of metastable 
vacua with n > 0. For vacua with k = Nf — N c , the only potential mode of instability is 
the tunneling of the n fourbranes through the potential barrier in the $ direction. The 
rate for this process can be suppressed by tuning g s and other parameters of the brane 
configuration. 

For k < Nf — N c , the brane configuration of figure 4 has another mode of instability. 
Some or all of the n D4-branes can connect with color D4-branes stretched between the 
NS and A^-branes and move towards the k D4-branes stretched between the NS and 
AS^-branes. Depending on whether n is larger or smaller than Nf — N c — k, this process 
leads to metastable vacua of the previous kind, or to supersymmetric vacua with a different 
value of k. 

The reconnection process described above is a tachyonic instability or a tunneling 
process, depending on the distance between the ends of the various D4-branes on the 
A,S{-brane in figure 4. That distance is in turn a function of the parameters of the brane 
configuration. It was shown in |j2j that in the geometric brane regime one can choose the 
parameters such that there is no tachyonic instability. In this regime, all vacua in figure 
4 are locally stable. Of course, the vacua with k < Nf — N c have shorter lifetimes than 
those with k = Nf — N c , due to the additional decay modes. All lifetimes can be made 
arbitrarily large by decreasing the string coupling g s . 

From the effective field theory point of view, the difference between the cases k = 
Nf — N c and k < Nf — N c is the absence or presence of the fields <p, <p in the decomposition 
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(j2.11|) . These fields correspond in figure 4 to strings stretched between the n and Nf—N c —k 
D4-branes. Their potential includes a tachyonic contribution at tree level, and a positive 
one due to the expectation value of the field $ n at the local minimum of its potential. In 
the brane regime, the potential for $ n is due to classical gravitational interactions, which 
are strong enough to push the system out of the tachyonic region [P2|,7| • In the field theory 
regime, where this potential is due to one loop effects in the low energy WZ model, the 



vacua with k < Nf — N c are always destabilized by (p condensation [14]. 



4. Early universe dynamics 

At early times, the brane systems described above are expected to be in an excited 
state. The question we would like to address is what happens to such a state as it evolves 
and, in particular, which of the vacua described in section 3 is approached by the system at 
late times. In this section we will address this question, first for the ISS system of figures 
1, 2, and then for the generalized ISS system of figures 3, 4. 

In a typical excited state of a brane system of the sort described in section 3, the excess 
energy can go into exciting the fivebranes and/or the .D-branes. It is natural to expect 
that for small string coupling, most of the energy goes into exciting the fivebranes, since 
they are the heaviest objects around. A generic excited state of an iVS'S-brane corresponds 
to a charged black brane. Thus, we will discuss the brane systems of section 3 with the 
NS and iVS'-branes taken to be non-extremal. i 

The non-extremality of the fivebranes leads to modifications of the potential energy 
landscape. Our main purpose below is to map out this landscape. Once this is understood, 
one needs to study the real time dynamics of the D4-branes in this landscape and, in 
particular, analyze the possible late time limits of the time evolution. We will briefly 
comment on this below, leaving a more complete study to future work. 

In the analysis []5|,|12"H of the brane systems of section 3, a dominant role was played 
by the gravitational attraction of the flavor L>4-branes to the iVS-brane. We will start by 
analyzing the modification of this attraction due to the non-extremality of the iVS'-brane, 
and then take into account some additional effects. 



5 DK thanks S. J. Rey for describing his unpublished work [23| on Hanany-Witten constructions 
with non-extremal fivebranes, in June 2008. 
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The background of an excited iVS'-brane, which wraps the IR 3 ' 1 labeled by (x\ t) and 
the v -plane, is given by the non-extremal fivebrane solution! |24j j: 

ds 2 = -f(r)dt 2 + H(r)[f(r)- l dr 2 + r 2 dQ 2 ] + (dx 1 ) 2 + dvdv , 

(4.1) 

Here (r, O3) are spherical coordinates on the IR 4 transverse to the iVS'-brane. The radial 
coordinate is given by r 2 = (y — i/ns) 2 + (x 7 ) 2 + ww. The functions /, H are given by 

f(r) = 1 - 4 > H(r) = 1 + 4 • ( 4 - 2 ) 

if is the dilaton field; e^ = g s is the string coupling far from the fivebranes. 

The geometry (|4.1|) describes a black brane with a horizon at r = r^. The size of the 
black brane (r^) increases with its energy density, £, as follows 

£ =wki( 1+ i)- (43) 

As rh — > 0, ( fO| ) approaches the extremal fivebrane (CHS) geometry [PS] , and ( |4.3|) ap- 
proaches the tension of a BPS N Sb-hiane. For ^> l s , the energy of the non-extremal 
fivebrane is much larger than its charge, and one can set H{r) ~ 1. This leads to the 
4 + 1 dimensional Schwarzschild black hole (extended in the five dimensions (x l ,v)). Its 
Hawking temperature is given as usual by Th = 1/r^. 

The non-extremal brane loses energy via Hawking radiation. As a result, its horizon 



shrinks at the rate 26 



j t (rl) - -9% ■ (4.4) 

If the initial radius is of order Z s or larger, eq. ( |4.4| ) implies that it will take a time t ~ l s / ' g 2 
or more to lose its excess energy. This is natural, since brane decay is a quantum process, 
and is thus suppressed by a closed string loop factor. The decay rate of a D-brane is 
proportional to g 2 ~ g s . The decay time of a D-brane is thus ~ l s /g s , which is much 
shorter than that for an iVS'S-brane (for j s < 1). 

The dynamics of the D-branes in figure 2 in the background of the extremal fivebrane 
( (|4.1| ) with rh = 0) was analyzed in 0. Our purpose here is to generalize this analysis to 
rh 7^ 0. We will split the discussion into two parts. In the next subsection we discuss the 
case [i = (see ( |2.1|) ); in the following one we move on to the system of figure 2, in which 



There is also a non-zero NS H field proportional to the volume form on the transverse S 3 . 
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4.1. n = o 

As discussed in section 3, the system with \i = in (|2.1| ) corresponds in the brane 
language to the configuration of figure 1, with the displacement of the iVS^-brane v<i set 
to zero. This is the magnetic SQCD brane configuration of |23] (with gauged flavor group 
0,0). It has a moduli space of supersymmetric vacua labeled by the expectation value of 
$ ( f2.1| ). In the brane construction, this moduli space corresponds to displacements in the 
MJ-plane of the Nf flavor .D4-branes 

Making the iVS'-brane non-extremal (i.e. turning on a non-zero (|4.2| )) breaks super- 
symmetry and generates a potential on moduli space. In order to compute it to leading 
order in g s , it is sufficient to consider a single L>4-brane stretched in y between the NS'- 
branes, and calculate its energy as a function of w in the background (fO|). 

The Dirac-Born-Infeld (DBI) action for a L>4-brane in the background of an iVS'-brane 
is given by 

/rV2 
d 4 x dye-*y/-detP(G) ab , (4.5) 

where T 4 j ' g s = 1/ {2Ti)H 5 s g s is the tension of the £>4-brane, and P{G) is the pullback of the 
spacetime metric ( |4.1| ) to the worldvolume of the brane. The potential is determined by 
taking w(x M ) to be constant, 
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V(w) = - dyJl- jj- - *' ' kJ9V> . (4.6) 



As in 0, we take the limit y^, Ay <C |yjvs|j in which the L>4-brane becomes a local probe 
of the geometry, and (|4.6|) simplifies to 



T 4 Ay / i r 2 h \w\ 2 T 4 Ay 





w 


2 


(vis + 1 


W 


2"J2 



9s V KVns + \ w \ ) 9, 



1 



1 rl 


w 


2 


2 (y 2 NS + 1 


W 


2)2 



(4.7) 



In the last step we expanded the potential to leading order in rh/yNS- We will restrict 
the discussion to the region l s <C <C Ay <C \vns\j an d work to leading order in the 
corresponding small parameters. This allows us to neglect some subleading effects. For 
example, we have taken the w-position of the D4-brane to be independent of y, i.e. assumed 
that the D4-brane continues to stretch in y when r/j > 0. One can show that corrections 
to this approximation are down by a power of rh/yNS- 

So far we assumed that the length of the D4-brane, Ay, is independent of w. This 
is certainly the case for = 0, in which the iV£"-branes are located at a fixed value of 
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y. However, in the non-extremal case they are deformed by the gravitational potential 
of the AS'-brane. To complete the calculation of the potential (|4.7| ) we need to take this 
deformation into account. 

The shape of an A£"-brane in the NS background (|4.1| ) is determined by a Nambu- 
Goto type action. The equations of motion of this action, which are analyzed in appendix 
A, lead to the shape 



V = V% + (Ui - Vns) 



4((y< - Vns) 2 + \w\ 



where yi gives the asymptotic position of the AS^-brane with i = 1, 2 in figures 1-4. The 
shape of the A£"-branes in the background of the AS'-brane makes the length of a D4- 
brane stretched between the two A£"-branes w dependent. Using (|4.8|) we find (to leading 
order in yi/yns) 



Ay(\w\) ~ Ay 



1 + 



Plugging this into (|4.7|) , we find the potential 



2 (Vns + \ w 
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(4.10) 



We see that the net effect of the non-extremality is to stabilize the L>4-brane at the origin. 
Some comments on the above analysis are in order: 

(1) We focused on a single D4-brane stretched between iV5"-branes in the non-extremal 
fivebrane background (4TT). Coming back to the brane configuration of figure 1 with 
t>2 = 0, the above analysis implies that for finite the Nf flavor branes are localized 
at w = 0. The supersymmetric moduli space is replaced by a vacuum with unbroken 
gauge symmetry, like in the corresponding field theory analysis. 

(2) We took the iV^'-branes, between which the D4-brane stretches, to be extremal. As 
mentioned above, it is more natural to take them to be non-extremal as well, with the 
same horizon size as the A^S'-brane. A iD4-brane stretching between two non-extremal 
AS"-branes with horizon sizes can be thought of as ending on the horizons. Thus, 
its length is smaller by 2rh than in the extremal case. This does not change the above 
discussion, since the resulting shift in energy is independent of w. 

(3) Although the non-extremality of the AS"-branes does not qualitatively affect the D4- 
branes, it nevertheless has an important effect on the dynamics of the system. For 
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finite r/j, the two parallel iVS"-branes attract each other and merge after a time that 
scales like (g s )°, much shorter than the scale of variation of (see ( |4.4p ). Thus, the 
fivebrane configuration of figures 1,2 is a very unnatural endpoint of the time evolution 
of this system.0 



4-2. ISS system 

We now turn to the case \i 7^ 0, which is described in terms of branes by the config- 
uration of figures 1, 2. Recall first the situation in the extremal case = 0. The brane 
system of figure 1 appears to have flat directions corresponding to the locations of the 
Nf flavor D4-branes in the u>-plane. However, gravitational attraction to the iVS'-brane 
provides a potential for these fourbranes, which localizes them at the origin |5|] . 

At that point the magnetic quarks, corresponding to strings connecting color and 
flavor branes, are tachyonic. Condensation of these tachyons corresponds to a reconnection 
process, which leads to the configuration of figure 2. The iV c flavor branes that remain 
between the iV^'-branes are still localized at the origin, due to the gravitational potential 
of the iVS'-brane. Thus, the resulting non-supersymmetric vacuum is locally stable. 

We would like to extend the above picture to the non-extremal case 7^ 0. Consider 
first the potential on pseudo-moduli space. Following |7j we take the limit Ay <C \yNs\, 
in which the shape of the L>4-brane in the (v,y) plane is well approximated by a straight 
line 

v = a(y - yi) . (4.11) 



We assume that the slope parameter, 



V2_ 

Ay 



(4.12) 



is small. The potential V(w) obtained from the DBI action (575) is given to leading order 
in r h /y NS and a by 
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(4.13) 



7 This conclusion can be avoided either if the iVS'-branes are initially at a very large distance, 
so that their motion towards each other takes longer than the time during which approaches 
zero, or if the rate of decrease of ru is in fact faster than that due to Hawking radiation. 
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Taking into account the bending of the iVS"-branes, (fOj), we find for \w\ -C \i/ns\ : 



V(w) 



T 4 Ay 



1+ v \ 4 " \wr + o(\ w \ 



(4.14) 



9a 

We see that, as one would expect, near the origin of pseudo-moduli space the potential is 
a sum of the attractive contributions due to the two sources of supersymmetry breaking 
- non-zero ( |4.2| ), and non-zero a, ([4.11 ). For ^> l s the contribution from non- 
extremality is much larger than the other one. 

The second element of the discussion of the extremal case that needs to be reexamined 
for rh 7^ is the reconnection of color and flavor branes, once the latter are localized at 
w = 0. If the iVS^-brane is extremal, the situation is as before; Nf — N c flavor branes 
connect to the color branes and slide off towards the configuration of figure 2. However, 
as mentioned above, early universe dynamics is expected to give rise to non-extremality 
of the A^S"-branes as well. Their horizon radii are likely to be comparable to that of the 
iVS'-brane; for simplicity we will take them to be equal below. 

Thus, we need to take the iVS^-brane to be non-extremal with horizon radius r^, and 
repeat the analysis of the tachyonic instability associated with the reconnection process 
mentioned above. This problem is studied in appendix B; here we summarize the results. 

For 

4a 3/2 Ay 4 

r h > r max ~ - = —=y/av 2 , (4.15) 



the configuration denoted by a long-dashed green line in figure 5, which consists of discon- 
nected color and flavor branes ending on the horizon of the iVS^-brane, is the only solution 
of the equations of motion for the D4-brane. In particular, in this range of rh, the brane 
reconnection instability discussed above is absent, and the gauge symmetry is unbroken. 

For rh < r max , two additional stationary points of the DBI action for the L>4-brane 
appear. These configurations, denoted by the solid red and short-dashed blue lines in figure 
5, correspond to a saddle point and local minimum of the energy, respectively. They break 
part of the gauge symmetry of the brane system. The saddle point (solid red line in figure 
5) has, by construction, a higher energy than the gauge symmetry preserving configuration 
(long-dashed green line), and local minimum (short-dashed blue line). When rh is in the 
range 

a? . a . A „. 

r max > r h > r rec ~ ~^^V = ^2 , (4.16) 

the symmetric configuration has lower energy than the symmetry breaking one. For rh < 
r rec , the minimum energy configuration is the symmetry breaking one. The energetics of 
this brane system is summarized in figure 6. 
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Fig. 5: Stationary points of the DBI action for a D4-brane for different value of 
ru- The short-dashed blue line corresponds to a gauge symmetry breaking local 
minimum, the solid red line is a saddle point, and the long-dashed green line a 
symmetry preserving local minimum. The parameters in the text were chosen as 
follows: a = 0.1, Ay = 0.l\y NS \, r h /v 2 = 0.2,0.136,0.025 for figures (a),(b),(c), 
respectively. 
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Fig. 6: Energies of the brane configurations of figure 5 as a function of Th (same 
parameters and color /dashing scheme as there). 



So far we discussed the effects of the non-extremality of the iVS , 5-branes on the ener- 
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getics of D4-branes stretched between them. In that discussion we implicitly assumed that 
the fivebranes remain static! as the D4-branes move in their background. As mentioned 
above, this assumption is generically violated. The gravitational attraction between the 
non-extremal iVS^-branes leads to collapse and merging of the fivebranes on a time scale 
much shorter than that associated with the change of . In that case, the endpoint of the 
dynamics is a system with coincident iV£"-branes, rather than that of figures 1,2. 

The analysis of energetics above is nevertheless useful, since much of it is applicable 
with minor modifications to the generalized ISS system, which will be discussed in the next 
subsection and does not suffer from the above instability. It may also be applicable to the 
ISS brane system of figures 1, 2 if the Schwarzschild radius decreases sufficiently fast 
to avoid the gravitational instability. We do not know why this should be the case, but 
for completeness briefly comment next on the early universe dynamics when the relative 
motion of the iVS^-branes can be neglected. 

Start at a time where rh is in the range ( |4.15| ). In this regime, the only minimum of 
the effective action for the £>4-branes stretched between A^5"-branes is the configuration 
of figure 1. Thus, the Nf — N c L>4-branes in figure 2 collapse to the configuration with 
unbroken gauge symmetry. 

As rh decreases with time, it enters the regime ( f4.16|) . A local minimum with broken 
gauge symmetry appears (short-dashed blue line in figure 5), separated from the symmetric 
vacuum by a potential barrier. In that regime, the broken vacuum has higher energy than 
the unbroken one (see figure 6). When further decreases below r rec , the symmetry 
breaking vacuum becomes the global minimum of the potential, but it is still separated 
from the symmetric one by a potential barrier. Since in the course of time evolution the 
system gets to this regime in the symmetric vacuum, the transition to the true vacuum 
can only happen via tunneling. The rate for this process is non-perturbative in g s . Thus, 
it does not occur on the time scale associated with the evaporation of the fivebranes, l s / ' g 2 s . 

As rh decreases further, it gets to the range rh < L, where the DBI analysis of 
.connection in appendix B bed. unre.iab.ef At that point, there are two plible 
scenarios for the time evolution. One is that the potential barrier between the two minima 

8 In field theory language, this is the assumption that the couplings in the Lagrangian can 
be taken to be constant in the early universe. In string theory it is natural for them to be 
time-dependent. 

9 Implicit in the above discussion is the assumption that r max ,^rec 3> h- Otherwise, the DBI 
analysis of the energetics becomes unreliable earlier. 
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in figure 6 remains finite for all > 0, and disappears in the limit. In that case, the 
tunneling rate increases with time, and eventually the system makes a first order phase 
transition to the symmetry breaking vacuum. The second possibility is that the barrier 
disappears for a finite value of r^, at which point the system makes a second order phase 
transition to the symmetry breaking vacuum. In both cases the endpoint of the dynamics 
is as in the field theory analysis ^-TL . 



Note that throughout the time evolution described above, the supersymmetric vacuum 
is separated from the ones studied here by a wide potential barrier. In the above discussion 
we implicitly assumed that the system starts at early times in the basin of attraction of 
the metastable vacuum. One can ask whether this assumption is natural. We will discuss 
this question in the next subsection in the context of the generalized ISS brane system, 
where the supersymmetric vacua are visible classically, and the gravitational instability 
associated with the separation of the iV£"-branes is absent. 

4-3. Generalized ISS system 

We start again by recalling what happens in the extremal case, = 0. The brane 
system of figures 3, 4 corresponds to the generalized ISS model (|2.9|), with the deformation 



parameter determined by the angle 9 via (|3.5|) . It has multiple vacua, labeled by the 
integers k, n in figure 4. Vacua with n = are supersymmetric; of the remaining ones, 
only those with k = Nf — N c are locally stable in the field theory regime (see section 2), 
and for general values of the geometric parameters in weakly coupled string theory (section 
3). 

As indicated in figure 4, for 9^0 the D4-branes stretched between iVS^-branes have 
a finite extent in w. Thus, one needs to be careful in parameterizing their locations in 
the w-plane. One way of doing that is to label the location of a fourbrane by the position 
of its endpoint on the iVS^ -brane, W\. The potential for such a L>4-brane can again be 
calculated from the DBI action. Expanding around the origin one finds 

n „ l) = Zi^f 1 _^ + ^>i + ... N ) . (4.17) 

9s \ w 2y% s J 



In (|4.17|) , Wo and w\ are taken to be real, although a priori they are arbitrary complex 
numbers. The location of the supersymmetric minimum, wq, can be chosen to be real and 
positive by a rotation in the (8, 9)-plane or, in field theory language, by a phase rotation 
on the chiral superfields. The potential for the complex scalar field w± has the property 
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that at the minimum its phase is the same as that of wq. Thus, for real wq, it is enough 
to restrict the potential for w\ to a real slice, as we do in ([4.171) . 
The potential ( 4.17|) has a minimum at w\ = u> m i n , 

«W„ * |F ■ (4.18) 

As 9 — > 0, wq — > oo, and the corresponding potential and minimum approach those of the 
ISS system. Since the potential ( 4.17 ) is valid only for w\ -C |j/jvs|, ( 14.18 ) is valid for 

„ \vns\ 3 ( . 1(V . 

W ° > 72 • ^ 4 " 19 ) 

This implies an upper bound on the deformation parameter 9, (pT5|) , ( p.6| ), 

0«P%- (4-20) 

\Vns\ 6 

We will restrict to this regime below, for simplicity. 

Another constraint on 9 that we will impose is associated with the requirement that 
all the vacua in figure 4 be locally stable. As discussed in |l^,0, the lightest modes of 
fundamental strings stretched between the color and flavor branes, which correspond in 
the field theory to the magnetic quarks denoted by ip, up in ( ^.11| ), are tachyonic if the 
distance between the endpoints of these strings on the AS^-brane, w±, is smaller than 
\/2nal s . Therefore, in order for all the vacua labeled by k, n, in figure 4 to be locally 
stable, one must have 

«* < ■ (4.21) 

The corresponding lower bound on the deformation parameter 9 (or ( p.5|) ) is easy to 
satisfy in the regime (|4.19| ), as is obvious from the fact that the corresponding bounds on 
^min, (namely y2iial s < «; m i n IvnsDj are compatible. 

Our main purpose here is to analyze the effect of non-extremality of the AS'S-branes 
on the energetics of the brane configuration of figure 4. As shown in appendix A, the 
introduction of non-zero 9 eliminates the instability of the ISS brane configuration noted 
above. Since the two AS"-branes are no longer parallel, the attraction between them does 
not lead to gravitational collapse. Instead, the fivebranes take a static shape in which 
each of them develops a localized dip towards the other (see figure 7). They also develop 
a logarithmic shape far from the intersection. Since we take the extra dimensions to be 
non-compact, this does not present particular difficulties. 

To understand the early universe dynamics of the brane system of figure 4, we next 
study the energetics of the D4-branes in the resulting fivebrane background. We discuss 
first the vicinity of the metastable minimum, and then that of the supersymmetric vacuum. 
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Fig. 7: Bending of iVS'-branes. Axes are offset, w — > w — wo, y — > y — yi — §Ay. 
This figure was generated using the results in appendix A, with yo = 1, Th = 0.1, 
and = 0.2. The expansion parameter rh/{6yo) = 0.5, which may be too large to 
be trustworthy, but was chosen to accentuate the features. 



For small w\ , a generalization of the calculation leading to (|4.17| ) gives 

T 4 Ay 
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(4.22) 



The minimum of (|4.22|) is located at 



Wr 



1 + r 2 /2a 2 l 2 



(4.23) 



with w m i n given by (|4.18| ). For small (ru <C al s ), w m i n (r^) is close to (|4.18| ). As 
increases, u> m i n decreases, and for al s it goes like w m i n (r^) ~ 2a 2 yj^ s /worf L . 

For w\ ~ wo, fc-e- in the vicinity of the supersymmetric vacuum, one expects on 
general grounds the non-extremality of the A^S"-branes to lead to important effects. The 
interaction of the fivebranes depends on two dimensionless parameters, 6 and r^/Ay. As 
discussed above, for 9 = and any finite r/j the branes collapse onto each other and form a 
single black brane with a larger horizon. On the other hand, for = they are mutually 
BPS and are not deformed at all. In general, one expects a smooth interpolation between 
the two limits. 
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This problem is discussed in appendix A. As increases, the branes bend towards 
each other in a small region near the intersection (see figure 7). The angle between them 
near the intersection also decreases as follows: 



= 



1 - , 2^) >+ 



(4.24) 



The expansion parameter on which the shape of the fivebranes depends turns out to be 
rh/OAy. When this parameter becomes of order one, the expansion breaks down. We 
expect the horizons of the two black branes to merge around that point, in a region in w 
whose size increases with ru- When that happens, any D4-branes stretching between the 
AS'-branes near wq disappear into the black hole. 

Note that the value of above which D4-branes near wq disappear from the dynamics, 
Th ~ 9 Ay, is very small. Indeed, using (|4.2U| ) one can show that 9 Ay <C al s . Thus, typically 



the supersymmetric minimum at w ~ wq is surrounded by a horizon throughout the time 
evolution of the system. 

We are now ready to discuss the early universe dynamics of the branes in the con- 
figuration of figures 3, 4. As before, we will take the Schwarzschild radii of the three 
A,S5-branes, r^, to be the same, and consider the dynamics as decreases. We start 
at an early time with ^> 9 Ay, so that any branes in the basin of attraction of the 
supersymmetric vacuum at W\ ~ fall into the black hole. Hence, all the branes that 
survive at late times must be located far from the supersymmetric vacuum, at a relatively 
small value of w. 

To follow the evolution of these branes, we start with in the range (|4.15[ ), so the 



D4-branes stretched between the NS' 2 and NS-branes are unstable to collapse towards the 
AS^-brane. Thus, at this point in the evolution, the system is attracted to a configuration 
in which the Nf — A c color D4-branes are stretched between the NS and NS[ -branes, and 
the Nf flavor branes are stretched between the A£"-branes at a value of w given by ( f4.23|) . 
As rh decreases, two things happen: 

(1) The position of the metastable minimum G4.23| ) increases. 

(2) The potential barrier that prevents the color and flavor D4-branes from reconnecting 
decreases (as discussed in the previous subsection). 

The two effects act in opposite directions - the first one tends to suppress the reconnection 
process, while the second enhances it. For ^> Z s the potential barrier prevents recon- 
nection from happening. What happens in the stringy regime is not clear, since the DBI 
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analysis of the reconnection process in appendix B breaks down there. What is clear is that 
one can choose the parameters such that the position of the metastable minimum ( |4.23| ) in 
the stringy regime is sufficiently large, so that reconnection does not occur as decreases, 
i.e. w m i n (rh) > \ / 2%al s for < l s . This is typically the case when the inequality ( |4.21| ) is 
satisfied by a wide margin. 

It might be that there is a small range of parameters for which w m [ n ( f4.18|) is suffi- 



ciently close to V27ra/ S that reconnection does occur for some < / s . Establishing the 
existence of such a range requires generalizing the analysis of appendix B beyond the DBI 
approximation; we will leave this for future work. 

Assuming that w m i n is taken to be large enough that the reconnection of color and 
flavor branes does not occur, the brane system approaches at late times the configuration 
of figure 4 with k = and n = Nf, in which the gauge symmetry is unbroken. This is one 
of the vacua that exists in the brane system but not in the field theory of section 2. 

If there is a region in parameter space in which the reconnection of color and flavor 
branes does occur, the system approaches at late times the configuration with k = Nf—N c , 
n = N c , the generalization of the ISS vacuum to non-zero We stress again, that if such 
a region in parameter space exists, it corresponds to a limited range of the parameters, for 
which the minimum w m i n is sufficiently close to the bound \/2-koI s . 

To recapitulate, we find that if the parameters of the generalized ISS brane system 
are chosen to be in the range in which all vacua in figure 4 are locally stable, early universe 
dynamics typically (i.e. for generic values of the parameters and initial conditions) leads 
the system to the most stringy and most symmetric of these vacua. One can fine-tune the 
parameters and / or initial conditions so that the system approaches one of the other vacua, 
but this is not generic. 



5. Discussion 

The work reported in this paper is part of a long term program of building phe- 
nomenologically viable models of nature in string theory, using configurations of inter- 
secting D-branes and A^iSS-branes. The first steps in this program were taken in the late 
1990's when it was shown that such configurations, which arise naturally in string theory, 
provide a simple way to realize supersymmetric gauge theories of the sort expected to play 
a role in nature. An added bonus of the string theory realizations of these theories was 
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an improved understanding of non-trivial dynamical phenomena such as mirror symmetry, 
Seiberg duality and Seiberg-Witten curves (see |1J for a review). 

In the last few years it became clear that these brane systems are also very useful for 
studying supersymmetry breaking. They exhibit metastable vacua of the sort discovered 
in the corresponding gauge theories in ||. Moreover, in some cases, such as the generalized 
ISS model, they have a much richer vacuum structure. The additional metastable vacua 
present in the brane systems but not in the corresponding gauge theories are in fact 
phenomenologically appealing for describing supersymmetry breaking in nature. 

In this paper we addressed the early universe dynamics of these brane systems. We 
found that they provide a picture which is closely related to the corresponding studies in 
field theory. The analog of turning on a finite temperature in the field theory analysis is 
replacing the A^S^-branes by non-extremal black branes with a finite Schwarzschild radius 
Th- The analog of the finite temperature second order phase transition associated with 
condensation of the magnetic quarks q, q is the transition associated with reconnection 
of the D4-branes (discussed around figure 5). The analog of the destabilization of the 
supersymmetric vacuum by thermal effects discussed in || is the merging of the non- 
extremal iV£"-branes at an angle 9 into a localized black hole that swallows any .D4-branes 
stretched between the iV5"-branes near the intersection of the two. 

The string theory realization provides again some unexpected bonuses: 

(1) R-symmetry breaking, which in the field theory discussions is important for phe- 
nomenological reasons, must be introduced in string theory for dynamical reasons. 
Indeed, the original ISS system which has unbroken R-symmetry, is an unnat- 
ural endpoint of early universe dynamics, due to the attraction of the iVS"-branes to 
each other .0 

(2) As mentioned above, the generalized ISS system, which has broken R-symmetry, has 
many metastable vacua that do not exist in field theory and are phenomenologically 
appealing. Interestingly, early universe dynamics drives the system to the most stringy 
of these vacua. This vacuum has the largest unbroken gauge symmetry, but at the 
same time it is the least stable of all the vacua labeled by (n, k) in figure 4. It is an 
interesting candidate for the supersymmetry breaking sector of a model of nature. 

A natural next step in the program described above is to use the brane models to realize 
inflation in string theory, and connect the inflationary period to the dynamics studied here 

10 The brane construction in which the iVS^-brane is replaced by Nf D6-branes suffers from 
the same problem. 
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(which we imagine to be post-inflationary). It would also be interesting to couple these 
models to a sector that realizes the MSSM in string theory and study the experimental 
consequences of the resulting models. For both purposes it is important to make the extra 
dimensions compact. 

Acknowledgements: We thank A. Giveon, M. Porrati, S. J. Rey and R. Wald for discus- 
sions. DK thanks the Weizmann Institute and the organizers of the Fifth Crete Regional 
Meeting in String Theory for hospitality during part of this work. JM would like to thank 
the organizers of the New Perspectives in String Theory workshop at the Galileo Galilei 
Institute for Theoretical Physics and Simons Workshop in Mathematics and Physics 2009 
for hospitality during part of this work. DK and OL are supported in part by DOE grant 
DE-FG02-90ER40560 and the National Science Foundation under Grant 0529954. JM is 
supported in part by a Ledley fellowship. AR is supported in part by a Department of 
Education GAANN fellowship P200A060226. 



Appendix A. Interaction between non-extremal fivebranes 

On several occasions in section 4, we encountered an A^S'5-brane in the geometry 
produced by another non-extremal fivebrane. In this appendix we will analyze the profile 
of an AS'S-brane caused by such a geometry; in particular, we will derive equations (|4.8| ) 



and (424). Strictly speaking, this analysis is valid for a probe extremal fivebrane in the 
background of a non-extremal one. However, one expects that the results are applicable 
to the case when both fivebranes are non-extremal, provided that their separation is much 
greater than their horizon radii. 

A.l. Linear analysis for a general angle 

We begin with the background of a non-extremal AS'S-brane, that wraps the IR 3 ' 1 
labeled by (x 1 , t) and the v -plane, and sits at y = x 7 = w = 0. This is equations ( pT~l] ) and 
(U) in the text0, which we rewrite here in a more convenient fashion: 

ds 2 = -fait 2 + H [{f~ l - 1) air 2 + aly 2 + (olx 7 ) 2 + alwolw] + (olx*) 2 + olvalv , 

(A- 1 ) 



11 Except that here we put the source at y = to simplify notation. 
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with r 2 = y 2 + (x 7 ) 2 + \w\ 2 and 



/(r) = l-^|, H(r) = l+ l l. (A.2) 



Let us place a probe fivebrane in the background ( |A.1|) . The probe brane wraps the 1R 3 ' 1 
labeled by (x l ,t), and when rn = 0, it follows a straight (complex) line: 

v = aw , y = yo , X7 = . (A. 3) 

For small the leading correction to this profile can be written as0 

v = aw + r 2 .v(w,w) , y = yo + r h y(w,w) , £7 = . (A. 4) 

In the following analysis we will restrict to the regime ^> l s so that we may set H = 1 
in O- 

The dynamics of the probe fivebrane is governed by the Nambu-Goto type action, 

S = -T NS J d 6 ae- 2 V-detP(G9 a6 , (A.5) 

where P(G) a b is the metric induced on the worldvolume of the brane. In general ( |A.5| ) 
should contain terms involving the Neveu-Schwarz S-field: P{G) should be generalized to 
P(G + B) in the Nambu-Goto term, and there should be a Wess-Zumino term coupling 
the potential Bq, where dB§ — *dB, to the worldvolume of the probe brane. The contri- 
bution of the .B-field to the Nambu-Goto term vanishes on the embedding ( |A.4|) , while the 
Wess-Zumino term is in general nonzero, but vanishes in the regime ^> l s we work in. 
Expanding ( |A.5|) in the powers of r^, we arrive at linear equations for y and v: 

y [(1 + a 2 )ww+ (a 2 - l)y 2 ] - A(yl + wwfd^y = , 
aw [(1 + a 2 )ww + (a 2 — l)yo] + 4(y(j + ww) 3 d w d^jv = . 

These equations have a unique solution which is regular at w = and has y(0) = v(0) = 0, 
and it leads to the following expressions for y and v: 

r l a2 ^fVo + \ w \ 2 





w\ 


2 


4yo(2/o + 


\w 


2 ) 



v = aw < 1 + 







■w 


2 


4 


w 


2 


L y 2 o + 


W 


2 



4yo V Vo 

+ (l + a 2 )ln 



2m„ Vo + H 



12 



Here we are looking for static configurations. 
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The properties of this solution at small and large values of a will be discussed in the 
next two subsections; here we make a few comments regarding qualitative features of the 
general solution. 

(1) Consider the small w behavior of the solution for y. Expanding around w = 0, we 
find 



y = yo 



2 1 \2 /I 14 

r h\W\ / \w 



(A.8) 



% 4 V yt 

When a < 1, corresponding to an asymptotic angle between the branes greater than 
45°, the probe brane bulges out near w = and bends around the non-extremal brane. 
When a > 1, so that the asymptotic angle is less than 45°, the probe brane dips down 
towards the non-extremal brane. (When a = 1, the next term in the expansion ( |A.8j ) 
shows that the probe brane dips down in this case as well). 

(2) While v approaches the straight line v = aw as \w\ — > oo, y on the other hand has a 
logarithmic running in \w\ (unless a = 0). However, this behavior is only noticeably 
different from an asymptotically constant profile on distance scales much larger than 
those we consider in this paper. 

(3) When both branes are non-extremal with equal horizon radii, each brane will of course 
have an equal and opposite effect on the other. In this case, ( |A.7| ) accounts for the 
net effect, to leading order in r^/yg. 

A. 2. NS' -branes in the background of NS: a = 

To describe the interaction between the A^S'-brane and AS"-branes, we set a = in 
the above analysis. Then the role of the probe brane is played by the AS"-brane, which 
sits at v = (or a constant), and stretches in w. Its profile in y is given by 



y = yo - A T\ ^ • (a.9) 





w\ 


2 


4y (yo + 


\w 


2 ) 



Moreover, since the configuration with a = has a symmetry v <-> —v, the branes should 
remain orthogonal to all orders in r^. 

The value a — is rather special: unlike the general result ( |A.7| ) , equation ( |A.9| ) does 



not have a logarithmic running at large values of w. Then, assuming the distance between 
branes in the y-direction, L, is measured at infinity, we find the relationlll 



yo - ^ , (A.io) 



13 



Here, and in the following, we drop higher order terms in r h /y and r h /L . 
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which leads to a modified version of (|A.9|) . 



1 + 



(A.11) 



4(L 2 + |«;| 2 ). ' 

To bring this expression into the form (|4.8|) , we recall that L = yi — yjqs (where i = 1, 2 
for NS[, NS' 2 respectively), and shift the origin of the y-direction so that the iVS'-brane 
sits at y = y NS : 



V = Vi + (yi- Vns)- 



(A.12) 



4((i/< - yNs) 2 + \w\ 2 ) 

A. 3. Intersection ofNS[ and NS 2 : large a 

To study the interaction of the A^5"-branes with each other, we first interchange v and 
w in (|A.1|) , so that the background describes the geometry produced by the A^-brane: 

ds 2 = -fdt 2 + H [(f- 1 -l)d? 2 + dy 2 + (dx 7 ) 2 + dvdv] + (dx 1 ) 2 + dwdw , 



e 2(^o) = H (fl 



(A.13) 



with r 2 = y 2 + (x 7 ) 2 + \v | 2 . If we interchange the role of v and w in (|A.3|) and ( |A.4j) as 
well, we can describe the NS 2 profile by identifying a = cot 9: 



w 



v cot 9 + r\w{v 1 v) , y = yo + r h,y( v , v ) 5 x 7 = . 



(A.14) 



Then it is clear that the general solution to the linearized equations of motion can be 
obtained from ( |A.7| ) by making the replacements v <-> w, a — > cot^: 



V = Vo 
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r h cot2e ln fyp + \ v 
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v cot 9 < 1 + 
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+ csc 9 In 



yo + H 
% 2 



(A.15) 



We are interested in very small values of 9, according to (|4.20|) . (|A.15| ) shows that y 
has a logarithmic behavior at large v, but we can define a separation between branes Ay 
by truncating this running at some large value of v, 



Ay = y 



+ 



-ci ~ ?/o + 



-ci 



(A.16) 



4y # 2 2/o " ^ 2 2/o 
where ci is some constant of order one. As the non-extremality parameter increases, 
while the separation between branes Ay is kept fixed, the minimal distance between branes 
yo decreases. Inverting ( |A.16| ), we find 



yo ~ Ay 



9 2 Ay 



ci . 



(A.17) 
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As the value of increases, the branes are also becoming more parallel. The expression 
( |A.15p for w shows that, while the angle between the branes approaches 9 at large values 
of v, at v = one finds a smaller angle 9q: 



1 



i + 



(A.18) 



Using (|A.17| ) in ( |A.18| ), we can reproduce ( [4.24| ). 

If the minimal distance between branes, yo, becomes of order r^, then the horizons 
of the NS[ and A^S^-branes touch, and these objects merge into a single black brane. 
Unfortunately the relation ( |A.16| ) breaks down long before this point. flA.16| ) follows from 
( |A.15|) , which is based on our linear analysis. In order for (|A.15|) , ( |A.16|) to be valid, we 
must require that rf l /(9y Q ) <C 1. Looking at the higher-order contributions to ( |A.16|) , one 
finds a general expansion 



Ay = y 



n=l 



0y 



In 



(A.19) 



Assuming that this expansion is invertible, 



2/o = Ay 



OO / 

n=l ^ 



8 Ay) 



In 



Ay g 



6Ay 



(A.20) 



we conclude that the branes connect if the function g takes a value rh/(Ay) <C 1. For a 
smooth function g this happens close to the point where g{x) vanishes, and this leads to 
the expression for the critical non-extremality parameter (r^) merge : 



merge 



x 6Ay 



g(x ) = . 



(A.21) 



For r h > (r h ) 

mergei the NS^- and A^^'branes merge in some region close to d — 0. 
To conclude this subsection, we observe that the critical value of given by (|A.2Q ) 
is equal to zero for parallel branes. This is consistent with the fact that such branes 
experience a finite force per unit length, so for any nonzero value of the non-extremality 
parameter they merge in finite time, and no static solution describing separated branes 
exists. 
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Appendix B. D4-brane dynamics in the non— extremal NS[ background 

In section 4.2 we discussed the dynamics of D4-branes in the geometry produced by 
a non-extremal NS[ , and found that the D-branes had an interesting vacuum structure if 
the non-extremality parameter, rn, is sufficiently large. In this appendix we will present 
the details of the analysis which led to this conclusion. 

The NS [-brane geometry is given by 

ds 2 = -f{7)dt 2 + (^(/(rT^? 2 + ^dhl) + (dx 1 ) 2 + dwdw , 

(B.l) 

e 2(^-^ ) = H ^ _ 

Here (r, O3) are spherical coordinates on the four-dimensional space transverse to the 
brane0; in particular, r = y/Jy — yi) 2 + (x 7 ) 2 + vv. The functions /, H are given by 

/(f) = l-0, H(7) = l+ l l. (B.2) 

Starting with the geometry ( B.l| ) , we consider a D4-brane which is stretched between 
the NS- and iVS^branes and extended in the directions (x l ). Assuming that the D4-brane 
is attached to the iVS'5-branes at points (v, w, y) = (0, 0, yNs) an d (v, w, y) = (v%, 0, 2/2)7 we 
conclude that the L>4-brane is localized at (w,lm(v)) = (0,0), and its profile is described 
by y = y(Re(v)). It is convenient to introduce polar coordinates (p, cf>) in the plane spanned 
by (Re(v),y), 

Re(v)+i(y-y 1 ) = pe i <i } , (B.3) 

and to parameterize the worldvolume of the D4-brane in terms of t) and the radial co- 
ordinate p. The profile of the brane in the Re(v)-y plane, <f>(p), can be found by minimizing 
the DBI action, 

S = -T A J d^xdpe^^-detPiG + B) a 

d 4 x J dp^/l + p 2 f(p)(d p <f>) 2 . 

The limits of integration in the p-coordinate will de determined below. The action ( B.4j 
leads to the equation of motion for 4>(p), 



9 



(B.4) 



P 2 f(p)d p 



<l + P 2 f{p){d p 



, (B.5) 



14 We introduce tildes in this expression to distinguish the spherical coordinates centered on 
the NS[ from those centered on the NS. The latter were used in the beginning of section 4. 



33 



and integrating this equation we determine the location of the D4-brane in a plane pa- 
rameterized by the polar coordinates (p, <fi) : 

MP) = ± T ; = ~ rldr = + ^ ■ (B.6) 

The integration constants, (po^o), specify the location of the point where the distance 
between D4- and iVS'J-branes reaches the minimal value, po- Although the integral entering 
( B.6| ) can be evaluated in terms of elliptic integrals F and K, 



MP) = ±V^— 1 {*T ( r l) - F farcsin ^, 4) } + </>o , (B.7) 
Po I VPo/ V P Po/J 

it is convenient to work directly with (|B.6| ). 

The integration constants, (po,(f>o), must be determined in terms of the boundary 
data, i.e. the location of the endpoints of the £>4-brane. In the setup depicted in figure 5, 
the .D4-brane is stretched between two iVS'S-branes, and its endpoints are located atHH 



7T / 7T V2 

(Pi,0i) = (Vi -Vns,-^) , (P2,02) = (^J(Ay) 2 + v%,- -arctan— ) . (B.8) 

One solution, satisfying the boundary conditions, consists of two disconnected components, 
each stretching radially outward from the horizon to the endpoints (pi,<fii) or (p2,<p2) 
respectively. This solution has d p (p = along the entire trajectory, corresponding to 
Po = r h in ( P-6| ). We refer to this configuration as the "disconnected" solution, and it is 
represented by the green long-dashed line in figure 5. 

To look for connected solutions with p > rh, we proceed as follows. Assuming that 
NS[ is near-extremal and the parameter a introduced in ( f4.11| ) is small, 

— <1, — <1, 7T-0 2 + 01<l, (B.9) 
Pi P2 

we can simplify the integral appearing in (|B.6|) . It is easy to see that the inequalities ( p.9| ) 
imply that <po is close to zero, and we begin by evaluating <p2 — </>o, which is close to ^. 
Taking the plus sign in (|B.6[) and shifting 0, we define 



M P ) = Mp) - 0o = / l( f\C h i P 2V • (B.io) 

V(P ~ r h)(P ~ Po) 



15 For now we ignore a possible sliding of the endpoint of the D4-brane that is attached to the 
iVS-brane. We will come back to this effect at the end of the appendix. 
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To analyze this expression, it is convenient to introduce dimensionless quantities £ and k: 



•p 2 -pl 



Po 



,2 ' 



K 



2 2 

Po~ r h 



(B.11) 



Rewriting equation ( [B.10|) in terms of these variables, we find 



(P) 



vWWTTT^) 



(B.12) 



Notice that (f>+(p) is an increasing function of p and a decreasing function of k. This 
function reaches its maximal value at p = oo and k = (which corresponds to po = °°) ? 



(oo) 



d£ TV 



(B.13) 



If the D4-brane begins at a finite p = p2, then there are upper bounds on p and po, which 
decrease the maximal angle §™ ax . To give ^>™ ax which is close to tt/2, both k and Pol Pi 
must be small; this requirement leads to the following inequalities: 



Th < PO < P2 



;b.i4) 



Expanding ( P-12|) for small re, we find 



<Mp) = arctan^(p) - 



ft 



£(p) 



1+^ 2 (P) 



+ arctan£(p) 



+ 0{k 2 ) 



In the approximation ( |B.14| ), (|B.15|) simplifies to 



+ («)=^+|-^-^+rtf4) + o/''- 

2 p 2 8 



Pi 



Similar manipulations lead to 

<f>-(pi) =<Po- 



_ po nrl 
L2 Pi 



Pi 



Pi 



(B.15) 



(B.16) 



(B.17) 



and combination of these two expressions gives the angle between the endpoints of the 
L>4-brane: 

Po + Po + tst\ 
Pi P2 4 pi' 



arctana = tx — (<p + — </>_) ~ 1 h A 



(B.18) 
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Equation (|B.18| ) can be used to determine po for the given values of (r^, pi, p%, a). The 
solution does not always exist, and when it does, it is not unique. First we notice that the 
right hand side of ( |B.18| ) is bounded from below: 



El + El + Y4>^ r 4, pI = ?-£wd- . (B.i9) 



The critical value po = p c corresponds to the minimal angl 

dmin — —. \ i (B.20) 

4 Pc 

for which ( |B.18| ) has a solution. Alternatively, one can fix a, p\ and p2] then ( JB.20| j can 
be interpreted as a relation for the maximal deviation from extremality which still allows 
for a connected solution: 

V 27tt J pi + p 2 

If a < a m i n (or > (r/ l ) 7na;r ) then there is no (connected) brane configuration satisfying 
the boundary conditions. 

For a > a m i n (or < (rh) max ) equation ( |B.18| ) has multiple solutions, and to analyze 
them, it is convenient to rewrite ( |B.18| ) in terms of p c and a m i n : 

i ^(ri)\UriV . (B . 22) 

3 a min \PcJ A \PcJ 

One of the solutions of this cubic equation has negative po and is not physical. Now 
we compare the energy of the two physical solutions with the energy of the disconnected 
straight branes. 

Starting with a general expression for the energy density which follows from the action 

(ED, 

£ = ^ / dpJl + p*f(d p <t))S , (B.23) 
9s J 



and substituting the solution (|B.6|) for 4>(p), we find the energy of the curved £>4-brane as 
a function of the minimal distance po : 



£ ^f\j^) 1 \^r d j ! Lz4) 1, \ (B , 4) 

9s J P0 VP - Po/ 9s J P0 VP - Po 1 
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Here and below we replace arctana by a, which is assumed to be small. 
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The straight-brane disconnected solution has po = and energy 

^straight = / «P H / <#> = [Pi ~ T h + p 2 ~ r hj 

9s Jr h 9s Jr h 9s 

Subtracting this expression from ( |B.24| ), we find 



(B.25) 



AS = £ - Sstraight = — \r h - p + I dp 

9 s I J Po 



Pi 



p ~ r h 
p 2 -pI 



2\ 1/2 



+ {pi -> P2} 



L _ po+ pj^A r /po == _^ 

9s \ h ° Po Ji Vu 2 - l(Vu 2 - 1 + y/u 2 - (r h /p ) 2 ) 



+ 



(B.26) 



+ {pi -> P2} ■ 



Eventually we want to look at A£ as a function of (pi, p 2 , rh, o), but we postpone the 
substitution of po = Po(Pi 5 P2, r /i, into (|B.26| ) until a later stage. 
The connected minimum is favorable when AS < 0; this implies 



1 + 



Po 



(£i ).(I2l) 



du 



v^ 2 ~^T(V^T+ ^u 2 - (r h /p ) 2 ) 



+ {pi -> P2} < 2 . (B.27) 



The left-hand side of this expression is a monotonically-increasing function of rh/po (while 
Pi/th and pil r h are kept fixed), so it can be equal to 2 for only one value of pq. The 
explicit construction will demonstrate that this value satisfies the inequalities 



r h < po < pi, p 2 



^B.28) 



Assuming that po and keeping only linear terms in Th/ Po, we can evaluate the 

integrals in the left-hand side of ( p.27| ): 



-a 


a/w 2 — 1 — li 


Pi/po 
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a/w 2 — 1 — W 


P2/P0 [ 






1 




1 J 



< 2 



(B.29) 



In the regime ( p.28| ), we can expand (|B.29|) in small parameters and find the condition 

(B.30) 



2 . 4pip 2 r /l _ 8 pi 
Po > 



pi + p 2 7T r h 

which justifies approximation po ^ r^. For future reference we also write an approximate 
expression for AS which is valid in the regime ( [B.28 ): 



P 2 r h \ 



AS = ^(po-r h )\^-^-^\ = ^(p -r h )\l-^ l 
9s { Po 2pi 2p 2 J g s po { 8 pi J 



;B.31) 
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We found that a connected configuration characterized by po is energetically preferred 
over the disconnected configuration as long as (|B.30|) is satisfied. It is convenient to rewrite 
this relation in terms of the geometric data (such as a, rh, pi, p-z) rather than po, which 
is determined dynamically. Substituting ( |B.30[ ) into ( |B.18| ), we find that reconnection is 
energetically preferred if a is larger than a critical value a rec , 



_ , rhjpi + P2) 4\/2 fj£ 

a > a rec = l\\ — - _ w — a min . [ti.sz 

P1P2 3V7r V r h 



For angles satisfying this inequality (so that a ^> a m i n ), equation ( |B.22|) has two solutions 
with positive po, 

3a / 8p c a 
Po - Pc- =Pcy > ( B - 33 ) 

I Q>min I dree ( Tfi\ „ .\ 

Po^Pc\ ^—=Pc\ ( 777T — <Pc, (B.34) 



3a V a \32 p 

and only the first solution obeys the condition (|B.30| ) characterizing a configuration with 



minimal energy. Thus, for given values of (a, pi, P2, r^) satisfying ( |B.32| ), equation ( [B.33| ) 
characterizes (through ( |B.6| )) the profile of a connected D4-brane which corresponds to a 
global minimum of energy. The profile with po given by (|B.34j ) describes a local maximum 
of energy which separates the global minimum from the configuration of straight discon- 
nected branes, which corresponds to a local minimum of the energy. As a decreases, the 



trajectories characterized by ( H.33j ) and ( H.34j ) become closer and their energies change 



in opposite directions (see figure 7). At a = a rec the minimum characterized by ( |B.33j ) 
and the disconnected brane minimum have the same energy. For a m i n < a < a rec , the 
connected minimum is a local minimum, while the disconnected configuration is the global 
minimum. Finally, when a decreases to a m i n , the maximum and local minimum degenerate 
to the same profile, with po = p c , and disappear. 

In the setting of section 4.2, one varies rh while keeping a fixed. To analyze this 
situation it is convenient to rewrite ( p.32| ) 



as 

,2 



rh < (r h )rec = ■/' 1/92a , < {r h ) max . (B.35) 
4(pi + P2) 

We conclude that the reconnected brane (with po given by ( B.33[) ) corresponds to the 
global minimum if < (r^) rec , to a metastable state if (r/ l ) rec < < {rh)max and it 
does not exist for r h > {rh)max- 
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In the preceding discussion the endpoints of the D4-brane were assumed to have fixed 
locations, ( P-8|) . While it is clear that sliding along the iVS^-brane increases the energy 
of the D4-brane, sliding in the Re (v) -direction along the iVS'-brane seems energetically 
advantageous. This is true, but it is easy to show that this effect leads to a subleading 
contribution to the energy. 

We conclude this appendix by rewriting expressions ( |rj.21| ) and QB.35 ) in terms of 
coordinates of the branes. To leading order in a and yi/yNS, we find 

V 27tt J y 2 - y N s 3V3n ^ B ^ 

, s _ a 2 (yi - yNs)(y2 - yi) _ a 2 Ay 

[ r h)rec — \ — . 

4(y2 - Vns) 4 
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